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Understanding HTS cuprates based on the phase string theory of doped 
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We present a self-consistent RVB theory which unifies the metallic (superconducting) phase with the half-filling 
antiferromagnetic (AF) phase. Two crucial factors in this theory include the RVB condensation which controls 
short-range AF spin correlations and the phase string effect introduced by hole hopping as a key doping effect. We 
discuss both the uniform and non-uniform mean-field solutions and show the unique features of the characteristic 
spin energy scale, superconducting transition temperature, and the phase diagram, which are all consistent with 
the experimental measurements of high-T c cuprates. 



The low-energy physics of high-T c cuprates 
is essentially described by a doped Mott insu- 
lator Q. The half- filling insulating phase can 
be well understood based on the Heisenberg 
model, where the bosonic resonating-valence- 
bond (RVB) state |!| gives a highly accurate de- 
scription of both short-range and long-range spin 
correlations. And the mean-field bosonic RVB 
theory || provides a very useful analytic frame- 
work over a wide range of temperature. 

The real challenge comes from the doped case: 
How the doped holes will influence the AF spin 
background and by doing so shape their own dy- 
namics. A simple generalization of the aforemen- 
tioned AF mean-field theory || to finite doping 
has indicated that the motion of doped holes gen- 
erally gets totally frustrated at the mean-field 
level such that an additional spiral spin twist 
must be induced in order to gain a finite kinetic 
energy. But this is only an artifact due to a mis- 
treatment of the singular phase string effect intro- 
duced by doped holes to be discussed below. We 
will show that the bosonic RVB ordering, which 
characterizes the short-range spin-spin correla- 
tions at half-filling, actually favors the hopping 
of doped holes and the RVB order parameter will 
still control the doped phase just like in the un- 
doped picture consistent with the original 

motivation for introducing the RVB concept by 
Anderson Q. 

PHASE STRING EFFECT. The t J 
model in the Schwinger-boson, slave-fermion rep- 



resentation has the form Ht—j = H t + Hj: 
H t = -tY J H lJ B Jl +H.c, 
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where Hi- 



fi fj > Bji 



and Af 



J2a biabj-cr. Here is a fermionic "holon" opera- 
tor and bi a is known as the Schwinger-boson oper- 
ator. At half-filling, only Hj exists and the mean- 
field state H is described by the bosonic RVB or- 
der parameter A s = (A*,-). Note that here in the 

slave-fermion decomposition, Cj ff = fjbi a (— <r)\ 
the Marshall sign is explicitly built in through 
(— a) 1 such that there is no sign a appearing in 
A|.. But Bji in the hopping term acquires a sign 
a, which would lead to (Bji) = in the RVB 
state, namely, the absence of the bare hopping 
for doped holes at the mean-field level. 

However, if one follows the hopping of a hole 
continuously without making the mean-field av- 
erage at first, then one finds Q that a sequence 
of signs (-1) x (+1) x (+1)... = (-1)^ will be 
picked up by the hole along its path c, where 
N% counts the total number of | spins being ex- 
changed with the given hole, originated from the 
simple sign a mentioned above. To illustrate the 
importance of such a phase string effect, let us 
take the one hole case as an example by express- 
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ing the total energy as follows 

Sk = i? o 7 -^E e " lk ' (r, " rj)M ^' ( 3 ) 

ij 

Mij = P[{c};Ni](-i) N t, 

{c},JVi 

with using the Brillouin-Wigner formula, where 
Eq denotes the ground-state energy with the hole 
being fixed at a given lattice site. The energy 
gain due to the hopping comes from the second 
term of (3). In My the summation runs over 
all the possible paths {c} of the hole connecting 
i and j. The significant fact is that the weight 
functional P can be shown to be positive-definite 
at -Ek < Eq . (The derivation is very similar 
to that of the single-electron propagator in the 
one-hole case [|J and the details will be presented 
elsewhere.) It implies that the phase string fac- 
tor (— 1)^ is the only source of phase frustra- 
tion which cannot be "repaired" since P > 0, 
and thus is expected to dominate the low-energy 
physics. The only additional phase effect when 
there are many holes will come from the exchange 
of any two holes due to the fermionic nature of 
holons in the slave-fermion representation. At 
small doping, the phase string effect obviously re- 
mains dominant. 

Thus the conventional slave-fermion scheme as 
a very successful framework at half-filling is no 
longer useful at finite doping if one does not know 
how to deal with the singular, nonlocal phase 
string effect [which basically counts how many 
i (or f by symmetry) spins exchanged with the 
hole during its propagation) hidden in the formal- 
ism. On the other hand, since the phase string 
effect induced by doped holes is always present in 
the ground state and can be "counted" in terms 
of the number of exchange between holons and 
spinons, it is possible to incorporate the singu- 
lar effect explicitly into the wavefunction like the 
fermion statistics. It turns out that one can re- 
alize this through a unitary transformation 
through which the t — J Hamiltonian is "reg- 
ulated" as the local singular sign a is "gauged 
away" . The resulting exact formalism is known 
as the phase string representation, where the RVB 



order parameter becomes 

A|,=^(e-^«)M,- CT) (4) 

(J 

while in the hopping term (01), Bji — > Bji = 
E ff (e i<M ?*) b)J> ia taaAHij = (e iA «) h\h r Note 
that the fermionic operator fi now is replaced 
by a bosonic holon operator hi after the trans- 
formation while bia — > biu which is still a bo- 
son operator. In this new "bosonization" for- 
malism, even though the singular sign a in the 
original Bji is "gauged away" , the topological 
effect of the phase string (the total phase for 
a closed path) is precisely tracked by lattice 
gauge fields A f tj and A^. Here A s v} = A s tj - ^ 

with J2c A tj = 1 / 2 J2iec (E CT <™1) fOT a closed 
path C and ^ D 0° = it per plaquette. And 

Ec4 = 1 / 2 Y,iec n t ( n t and n t are s P inon 
and holon number operators, respectively.) Here 

Afj and A\j describe quantized flux tubes bound 
to spinons and holons, respectively. And </>y de- 
scribes a uniform flux with a strength tt per pla- 
quette. 

MEAN FIELD THEORY. In the phase 
string formulation, the hopping term becomes fi- 
nite as (6^) ^ H when A s ^ 0. Namely, the 
single RVB field A s will still be the only order 
parameter that controls both the undoped and 
doped phases. Note that A s decides the short- 
range spin-spin correlations as shown in the rela- 
tion (Si-Sj) = — I/2|A S | 2 for nearest-neighboring 
i and j. 

In the following let us focus on the so-called 
uniform-phase solution where holons and spinons 
are assumed to distribute homogeneously with- 
out the phase separation. In the ground state, 
the bosonic holons will experience the Bose con- 
densation (BC) such that one may replace A^ 
by A^j satisfying ^ D A^ = Sir (S is the doping 

concentration) and A{- — > —<j>% as Afj diminishes 
due to the pairing of spinons with a gap opening 
(see below). Then the mean-field state can be 
straightforwardly obtained. 

Fig. 1 illustrates how the spin dynamics evolves 
into the metallic regime: the spin dynamic sus- 
ceptibility at (7r,7r), x"7r( w ); i s shown at S = 
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and 6 = 1/7 w 0.143, respectively. Compared to 
half-filling, a resonance-like single peak emerges 
at a finite energy E g with the weight diminishes 
below it. The inset shows the doping-dependence 
of Eg which peaks at ~ 0.4J (independent of t) 
around 5 ~ 0.2, which is very close to the famous 
41meV peak observed || in optimal YBCO com- 
pound by neutron scattering (for J ~ lOOmeV). 
Eg decreases monotonically with reducing S in 
low-doping regime. 

The phase string effect is solely responsible for 
such a resonance-like peak through the gauge field 
A^j. Another unique prediction of the phase 
string effect is shown in Fig. 2 for the local 
susceptibility Xl( w ) — l/A r X} q x"( t l) w ) which 
shows multi-peak structure at higher energies like 
Landau levels as A% describes uniform fictitious 
flux seen by spinons. Note that the high-energy 
peaks are contributed by those spin fluctuations 
away from the AF momentum (7r,7r) in momen- 
tum space and are only significant under the inte- 
gration over the whole Brillouin zone, which re- 
mains to be seen experimentally in the uniform 
phase like the optimal-F-BCO compound. 

Another interesting property is that the ground 
state is intrinsically a superconducting state with 
d-wave symmetry. This can be verified by the 
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Figure 1. Dynamic spin susceptibility at the AF 
vector (7r,7r) as the function of energy. 
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Figure 2. Local dynamic spin susceptibility as the 
function of energy 



paring order parameter || 



(5) 



for nearest-neighboring sites i and j. Here the 
numerical factor changes sign from j — i ± x 
to j = i ± y determining the d-wave symmetry, 
ho = (hi), and the phase $| is defined as 



(6) 



$^ = Im In (zi 

l^i a 

which describes vortices (anti-vortices) centered 
at up (down) spinons. The vortices-antivortices 
described by $| are all paired up in the ground 
state to ensure a finite Af„- C . But at finite tem- 
perature, free vortices appear in <frf as spinons 
are thermally excited will eventually kill the su- 
perconducting condensation. Using the effective 
Hamiltonian for holon = —th "Yluj\ Hij +h.c, 
T c can be estimated [|J as the temperature at 
which the excited spinon number equals to the 
total holon number. Fig. 3 show the correla- 
tion between T c and the spin characteristic en- 
ergy Eg which is consistent with the experimen- 
tal relation for the cuprate superconductors |t],|| . 
In particular, T c ~ 100 A' at S = 0.143 with 
Eg ~ 40meV without tuning parameters except 
for J ~ lOOmeF. 

The phase diagram of the present mean-field 
theory is illustrated in Fig. 4. Note that the 
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whole theory is underpinned by the RVB order 
parameter A s controlling the nearest-neighboring 
AF spin correlations which are gradually reduced 
with the increase of S. Both AF long-range order- 
ing in the insulating phase (holes must be local- 
ized) and superconducting condensation in metal- 
lic phase occur inside such an RVB background, 
characterized by the BCs of spinons and holons, 
respectively. Note that there is an underdoped 
metallic regime denoted by SBC in Fig. 4 un- 
der the dashed curve which indicates the case of 
spinon BC in the metallic phase, corresponding 
to a non-uniform solution of the present mean- 
field theory ||. In the above uniform phase, the 
spin gap Eg shown in Fig. 1 prevents the BC of 
spinons. But a non-uniform distribution of holons 
will lead to a spatial redistribution of the ficti- 
tious flux seen by spinons through to gen- 
erate a low-energy tail in spinon spectrum. At 
low doping, the Bose condensation of spinons can 
happen as the tail extends to zero, which phys- 
ically means short-range AF orders exist in the 
hole-dilute regions. Whether such a solution is a 
micro or macro phase separation still needs fur- 
ther investigation. But generally the dynamic 
spin susceptibility function will show a double- 
peak structure in energy space near (tt, tt) with a 
pesudo-gap feature as there always will be some 
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Figure 3. Relation between the spin energy E g 
(Fig.l) and T c . 




Figure 4. Phase diagram of the RVB state char- 
acterized by A 3 : AF and SC phases correspond 
to spinon and holon Bose condensations, respec- 
tively. And SBC denotes the spinon Bose con- 
densation in the metallic regime. 



weight extends towards the zero energy limit 
This fact may have already seen experimentally 
in the underdoped YBCO compounds [@,||. 

REFERENCES 

1. P. W. Anderson, Science 235, 1196 (1987). 

2. S. Liang, B. Doucot, and P. W. Anderson, 
Phys. Rev. Lett. 61, 365 (1988). 

3. D.P. Arovas and A. Auerbach, Phys. Rev. 
B38, 316 (1988). 

4. Z. Y. Weng et ai, Phys. Rev. B55, 3894 

(1997) ; D. N. Sheng et al, Phys. Rev. Lett. 
77, 5102 (1996). 

5. Z. Y. Weng et al, Phys. Rev. Lett. 80, 5401 

(1998) ; Phys. Rev. B59, 8943 (1999). 

6. H. F. Fong et al, Phys. Rev. Lett., 75, 316 
(1995). 

7. P. Bourges, in The gap symmetry and fluctu- 
ations in high termperature suerpconductors, 
edited by J. Bok et at. (Plenum Press, 1998). 

8. H. F. Fong et ai, 3ond-ma t/9910041| . 



